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- ■ ■ Abstract. A cluster expansion is proposed, that applies to both continuous and 

, discrete systems. The assumption for its convergence involves an extension of the 

' neat Kotecky-Preiss criterion. Expressions and estimates for correlation functions 

^SJ , are also presented. The results are applied to systems of interacting classical and 

quantum particles, and to a lattice polymer model. 
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^ i 1. Introduction 

• Cluster expansions were introduced at the dawn of Statistical Mechanics for the study 

! of high temperature gases of interacting particles. They constitute a powerful perturbative 

' method that is suitable for geometrically large systems, such as encountered in Statistical 

^ ■ Physics. Numerous articles have contributed to the subject; a list of relevant publications 

I is |GKl ICaml Bry[ IKPT IPhl IDobI IBZ[ IMirl \baic\ and references therein. Cluster expansions 



O ■ can now be found in standard books, see Chapter 4 of Ruelle |R.uej . or Chapter V of 

Simon |Simj . They apply to continuous systems such as classical or quantum models of 
interacting particles, and also to discrete systems such as polymer models, spin models 
r~| ■ (with discrete or continuous spin spaces), or lattice particle models. The methods for 

"j^ i treating these various situations share many similarities, but a somewhat different cluster 

expansion was so far required in each case. 

The exposition of the cluster expansion is often intricate, and the paper of Kotecky and 
Preiss |KPj should be singled out for proposing a clear and concise theorem, that involves 
^ ' a neat criterion for the convergence of the expansion. This theorem applies to discrete 

^ . systems only, and its rather difficult proof was subsequently simplified in |Dob| IBZ| IMirj . 

Explicit expressions for the contribution of cluster terms can be found in the article of 
Pfister |Phj ; this helps clarifying the situation and allows for computations of lowest order 
terms. The condition for the convergence is the same as in |KPj in the case where polymers 
are subsets of a lattice. 

The goal of this paper is to present a general theorem that applies to both continuous 
and discrete systems. The condition for the convergence is given by an extension of the 
criterion of Kotecky and Preiss, see Equation © below, and the contribution of cluster 
terms involves explicit expressions. Furthermore, we derive expressions and estimates for 
correlation functions. 

The theorems presented here are illustrated in Section ^ in three different situations, 
namely classical and quantum gases of interacting particles, and lattice polymer models. 

2. Cluster expansions 

Let {A,A,fi) be a measure space; ^ is a complex measure and |/x|(A) < oo, where |^| 
is the total variation (absolute value) of Let C be a complex measurable symmetric 



2 DANIEL UELTSCHI 

function on A x A. The partition function Z is defined by 

Z=Y.^jMAi)...[df,{A^) n {^ + C{A,A,)). (1) 

n > 1 ^ i<i < n 

The term n = of the sum is understood to be 1. 

We denote by Gn the set of all (unoriented) graphs with n vertices, and Cn C Qn the set 
of connected graphs of n vertices. We introduce the following combinatorial function on 
finite sequences (^i, . . . , An) of A: 

, I 1 if n = 1 

V.(^i,...,A„)= (2) 
Z^GeCn ll(i,i)eG'^vA'^jJ n n ^ 2. 

The product is over edges of G. A sequence {Ai, . . . , A^) is a cluster if the graph with n 
vertices and an edge between i and j whenever (^{Ai,Aj) ^ 0, is connected. 

The cluster expansion allows to express the logarithm of the partition function as a sum 
(or an integral) over clusters. 

Theorem 1 (Cluster expansion). 

Assume that |1 + C(^,^')| ^ 1 for all A, A' € A, and that there exists a nonnegative 
function a on A such that for all A £ A, 

J dMA')\C{A,A')\e<^'^ ^a{A), (3) 

and J d|^|(A) e'^^'^) < oo. Then we have 

Z = exp| / dn{Ai)... / dfi{An)^p{Au...,An)y 

Combined sum and integrals converge absolutely. Furthermore, we have for all G A 

l+Y^nj d\^l\{A2)... I dl/il(A„)|(^(Ai,...,A„)| ^ (4) 

The rest of the section is devoted to the proof of this theorem; the reader interested in 
results only should jump to Section IHl that discusses correlation functions. 

Proof. Another inequality turns out to be helpful. Multiplying both sides of Q by 
|C(^, ^i)| and integrating over Ai, we find using 



j d\^i\{A^)... I d\^Ji\{An){J2\aAAi)\)\^{Au...,An)\^a{A) 

n>l-^ 1=1 



(5) 



for all Ae A. 

The strategy is to show inductively that (jSJ implies (jlj. Convergence of the cluster 
expansion follows and allows to prove Theorem ^ 
We prove that the following holds for all N, 

N 

1 + / d|/i|(^2)... / d\fi\{A„)\ipiAi,...,Ar,)\ ^ e"(^i). (6) 

n=2 •' •' 

The case = 1 is clear and we consider now any N . The left side is equal to 

N 

l + E/dl/^I(A2).../dl/.I(A„)^-^|^ W aA^,A,)\. (7) 
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Let us focus on the sum over connected graphs G. Removing all edges of G with one 
endpoint on 1 yields a possibly disconnected graph G' . Let (Gi, . . . ,Gk) be a sequence 
of connected graphs where Gi has set of vertices Vi, Vi U • • • U = {2,... ,n}, and 
Vi nVj = ^ ii i ^ j . Each sequence determines a graph G' , and to each G' corresponds k\ 
such sequences. See Fig. ^ for an illustration. Therefore 




Figure 1. Illustration for G, G', and (Gi, . . . , Gk). 



|E n c(^^'^.)|^ e n{ n c(a.,a,)2: n 

GeC„(i,j)eG k^i (Gi,...,Gfc)£=i (i,j)eG, G^ (i,i)eG^ 

(8) 

The sum over G'^ runs over nonempty sets of edges with one endpoint on 1 and one 
endpoint in (i = 1 in the last product). We have 



n aA,A,)=ll{l + CiAuA,))-l. 



(9) 



Using 



n-l 



J](l + Oi) - 1 = [J] (1 + a.) - l] (1 + an) + an, (10) 

i=l 1=1 

one easily proves inductively that the absolute value of © is smaller than X^jgy^ lC(^i) 

The sum over sequences (Gi, . . . , G^) can be done by first choosing the respective num- 
bers of vertices nii, . . . ,mk whose sum is n — 1, then by summing over partitions of 
{2, . . . , n} in sets Vi, . . . ,Vk with \Vi\ = mj, and finally by choosing connected graphs for 



each set of vertices. The number of partitions is 



mi\...mk'- 



Then ((Tj) can be bounded by 



1 + 



n=2fc^l ■ mi,...,mfc>l £=1 

me 

HA[,...,A'^^)\Y^\CiAuA[] 



irii-\ \-mk=n~l 



i=l 
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We can sum over n; the constraint mi + • • • + ^ — 1 can be relaxed into mi ^ N — 1 
for all i. Using © with n ^ — 1, we obtain the bound 

1+ j;i[a(^ir=e'^(^^). (11) 



This proves inequality @ . Absolute convergence of the cluster expansion follows from 
and summability of e"'-"^-' . 

The rest of the proof is standard and consists in showing that clusters are indeed the 
terms of the expansion of the logarithm of the partition function. The idea is to expand 
Z so as to recognize an exponential. 

Z = l+Y.k f^t'^^^y-- [MAn)Yl n ^(^-^i)- (12) 

The sum over n converges absolutely. We proceed as above and consider sequences of con- 
nected graphs (Gi, . . . ,Gk) whose sets of vertices form a partition of {1, . . . , n}. Summing 
first over the number of vertices of each partition, then over partitions, we get 

1 



k\ ^-^ mil . . . mfc! 

n ^ 1 fc ^ 1 mi,...,mj. ^ 1 
mi -I \-mk=n 



k 

n{ / MAi) ... [ dfi{Am,) Yi n c{A,A,)} (13) 

k 

E fljdfi{Ai)...jdfi{A^,)ip{Au...,Am,). 



n ^ 1 fc ^ 1 mi,...,mf; ^ 1 £=1 ' 
mi -I |-mfe=n 



Absolute convergence of the clusters allows to remove the sum over n, and this completes 
the proof. □ 

3. Correlation functions 

An advantage of the cluster expansion is to characterize correlation functions. The 
relevant general expressions are 

Z{Au...,Am) = (n-m)l f M^m+i) . . . ! M^n) n (1 + C(^.,^i)); 

(14) 

Z{Ai,...,Am)= V , / d^(A™+i)... / d^{An)^{Ai,...,An). (15) 

^-^ in — my. I I 

It is understood that in both expressions, the case n = m corresponds to taking the 
integrand without integrating on polymers. 

Theorem 2 (Correlation functions). 

Under the same assumptions as in Theorem^ we have 

k 



{Vi,...,Vfe}i=l 
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where the sum is over partitions of {1, ... , m}, i.e. ViU- ■ -DV^ = {1, . . . , m}, and ViCiVj = 

ifiy^j. 

The next result deals with estimates of correlations. To exhibit a suitable decay, an 
efficient strategy is to establish the criterion Q in a stronger form. We consider a non- 
negative function b on A, and a nonnegative symmetric function c on A x A (both can be 
identically zero, but the larger they are the better). We introduce the notation 

UA,A') = C{A,A')e<^''''^ 

c{Ai, . . . , An) = mill c{Ai, Aj) if n ^ 3 

if,{Ai, ...,An) = ^{Ai, ...,An) e^(^i'-'^") . 

The utility of functions b and c will be illustrated in Section ^ The following theorem 
contains estimates on correlations; compare with the definition (|15|) of Z(Ai, . . . ,Am)- 

Theorem 3 (Decay of correlations). 

Assume that |1 + (^c(^)^')l ^ 1 /''^ A, A' € A, and that there exists a nonnegative 
function a on A such that 

''d\fi,\{A')\UA,A')\e<^"^ ^a{A) (16) 

for all A ^ A. Then the following estimate holds true for all m'^ \, and all Ai, . . . , Am G 

A, 

V , [ d\fib\{Am+i) . ■ ■ [ d\nb\{An)\ipciAi,...,An)\ 

^ (n-m)\J J 

m 

^ exp{^[(l + 7)'"-l]5^a(^.)} n {l + \Cc{A,A,)\). 

i=l 1 ^ i<j ^ m 

Here, we set 

n 

7 = sup sup \\\{l + aAo,Ai))-l (17) 



i=l 



(clearly, ^ 7 ^ 2^. 

Notice that the case m = 1 is 



l+J^nj d\fib\{A2) . . . j d\^lb\iAn)\MAl,■■■,An)\ ^ e'^(^i) (18) 

which is reminiscent of @. Multiplying both sides of (fTH|l by |(^c(^)^i)| and integrating 
over Ai, we obtain from (|16|1 



f d\fib\iAl) . . . /"d[/ib[(A„)(^[Cc(A^i)|)!</'c(^l,...,^n 

n > 1 i=l 



)\^a{A) (19) 



for all A G A. H18() and ()19() can be proved exactly the same way as @ and ©. 

We turn to the proofs of Theorems and 01 Readers interested in applications should 
jump to Section 13 We start with the proof of Theorem |31 as it will imply the convergence 
of the terms appearing in Theorem |2J We assume that H19|) has been established. 
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Proof of Theorem\^ A connected graph G on {1, . . . , n} can be decomposed into a graph 
G' on {1, ... , m}; a partition {Wi, . . . , Wk} of {m+1, . . . , n}; connected graphs Gi, . . . ,Gk 
where Gi has set of vertices W^; non-empty subsets Vi, . . . ,Vk of {1, . . . , m}; non-empty 
sets of edges between Wj and each vertex of Vj . This is illustrated in Fig. |2 Conversely, 
choosing these graphs and sets of vertices yields a graph on {1, . . . , n}; it is not necessarily 
connected, but we obtain an upper bound by dropping this constraint. Therefore we get 




Figure 2. 



ggc„ {i,j)eG 



^c{Ai,...,A„) 



CGSm (i,i)eG' 



{Wi,...,Wk} 



n I E n 

£=1 ^ Ge&CiWe) (i,j)eGi 



exp{ mm c{A,,A,)}l[\ll{l + CiA,,A,))-l]. (20) 
y£C{l,...,m} ieVe j€We ^ 



From the definition (|17jl of 7, we have 



n|n(i+c(^-^^-))-i 



^ WE n(i+c(^-^.-))-i 



1^1 



We used © and ()1U|) for the second inequality. Furthermore 

E ^ E ic(^^^^.)i = -[(i+7r-i]EEK(^^'^.)i. 



(21) 



(22) 
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Then 

^ E n E Em E 

G'eSm {i,j)eG' n^mfc^O ni,...,nfc > 1 

riiH \-nii=n—m 



^c{A[,...,A'„^) 



Yi\^ [ diM.K^'i).-- / di/.6iKj| E n 

-. m ni >. 

— [(i+7r-i]EEiCc(^-^.)i 
1=1 j=i 



(23) 



The case A; = of the right side should be understood as the case n = m of the left side. 
The expression inside brackets in the right side converges because of p9j) : this allows to 
sum over n, and we get Theorem IHI □ 

Proof of Theorem\^ We expand the product in (|14() and we obtain a sum over graphs of 
Gn- A graph of Qn can be written as a sequence of graphs (Gi, . . . ,Gk,G'), where each 
Gi is connected and has at least one vertex in {1, . . . ,m}, and the vertices of G' are in 
{m + 1, . . . ,n}. 

Let Wi be the set of vertices of Gi, Vi = Wi H {1, . . . , m}, and W be the set of vertices 
of G'. {Vi, . . . , Vfc} is a partition of {1, ... , m} and W C {m + 1, . . . , n}. Furthermore, 
let nj = \Wj \ Vj\, 1 ^ j ^ fc, and p = \W'\. 

The number of partitions {Wi, . . . , Wk, W'} of {1, . . . , n} corresponding to given ni, 
Uk, p, and {Vi, . . . , V^}, is equal to ^fni^' Therefore 



{Vi,...,V^fc}n m ^0,...,nfe ^0 £=1 J J 

"iH hfife ^ n—m 

V.((A,)^ey,,^'i,---,0]^/dM^'i)---/dM^p)E n ^(^^'40' (24) 

GeSp (i,j)eG 

where p = n — m — X]i=i '^i- terms converge absolutely and uniformly in n because of 
Theorems n and El this allows to perform the sum over n, and we obtain Theorem |21 □ 



4. Illustrations 

The first two examples are classical interacting particles and lattice polymers — they 
are well-known but nevertheless constitute nice illustrations of the results above. The last 
example deals with quantum interacting particles and is more involved. 

4.1. Classical interacting gas. We consider a gas of particles that are subject to pair 
interactions. Let A = A be a bounded subset of M'^ and d^{x) = zdx, where z G M+ 
is the fugacity and dx is the Lebesgue measure. Let /3 be the inverse temperature and 
U{x — y) ^ represent the interactions between particles at positions x,y € A. Setting 
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C{x,y) = e ^^^^ ^) — 1, the partition function is given by (^. We choose a{x) = 1 and 
we easily check that the criterion Q holds whenever 

z [ dx{l - e-^^(^) ) ^ e-i . (25) 



This condition for the convergence of the cluster expansion is well-known, see |Ruej Chap- 
ter 4. One then obtains an expression for the thermodynamic pressure p{(5,z), namely 

/3p(/3, z)=z+Y, ^"^^ / "^^^ • • • / ^1' • • • ' (26) 



n > 1 



where integrals are over M^. This expression is absolutely convergent because of Q. 
Furthermore the function ip{-) is analytic in /3; by Vitali convergence theorem (see e.g. 
|Simj Theorem V.2.7 in the context of statistical physics), p(/3, z) is analytic in /?, z when 
the condition (^5]) is satisfied. 

We study now the correlations. For bounded A, we consider the functions 



■y ^ z"' f f -I— r 

Plixi) = ^ ^ -r—^ dx2 . . . dxn Yl {^ + C{xi,xj)), 

2 ^ z"' f I s—r 

P2{xi,X2) = -^-r^-—r^ dx3... dxn n {l + CiXz,Xj)). 
By Theorem|21 the truncated two-point correlation function p\[xi,X2) is given by 

^2(3:1,3:2) = P2{XI,X2) - Pl{xi)pi{x2) 

= ^ n{n - l)z" dx3... dx„ (f{xi, x„). (27) 

n ^ 2 

We use Theorem 131 (with 7 = 1) to get a bound for the decay of correlations. Let 5 = 
and c(x) ^ satisfy the triangle inequality, such that 

z [ dx{l - e-^^(^) ) e'^^^') ^ e-i . (28) 



We have 

e'^^'-'-'^\pl{xi,X2)\ ^ Yl ^{ri-l)z^ I dx3... / dx^^i, ■ ■ ■ ,Xn)\e<^'-^-^ 



n>2 



^ e^ {1 + (1 - e-PV{xi-x2) ^ gc{xi-x2) |_ (29) 

The right side converges to e^ as \xi — X2I ^ 00. This shows that 

\p\{xi,X2)\ ^ const 6-^(^1-^2) (30) 
for all functions c satisfying (|28|) . 



4.2. Polymer models. A polymer is a connected subset of 7/. Let A be the set of 

polymers in a finite set A C 'L'^. The measure p is taken to be the counting measure 
multiplied by a weight w{A) satisfying |u)(A)| ^ e"^'"^' with rf = 2\og{2d(f)) + Here 

(j) = is the Golden Ratio. Polymers interact through a condition of non-intersection, 

that is, C(^, ^') is -1 if vl n ^' / 0, and is otherwise. 

To check the criterion Q, we choose a{A) = (j)~^\A\. It is enough to consider the case 
where A = {0}. If ^ is a connected set, there exists a closed walk with nearest-neighbor 
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jumps whose support is A, and whose length is at most 2\A\. This can be seen by induction: 
knowing the walk for A, it is easy to construct one for A U {x}. The number of connected 
sets of cardinality n that contain the origin is therefore smaher than the number of walks 
of length 2n starting at the origin, which is equal to (2d)^". The left side of (jSJ is bounded 
by Y.n ^ e-iri-<t>-^)n and this is equal to (f)"'^. 

The cluster expansion provides absolutely convergent series for thermodynamic quanti- 
ties. Many physical models can be mapped onto a polymer model, and TheoremOltypically 
provides informations on correlation functions. 

4.3. Quantum interacting gas. The description of a gas of quantum particles in a 
bounded domain A C M"^ should start with the state space of the system, which is the Fock 
space !F± = 0^ -> q P-j-(L^(A))'^^, where P+ (resp. P_) is the projector onto symmetric 
(resp. antisymmetric) functions. Next one should introduce the Hamiltonian (Laplacian 
for the kinetic energy of the particles, and another operator for the interactions). One 
could then write down the partition function Z = Tr . 

It serves better our purpose to define the model in the Feynman-Kac representation. 
See |Ginj for a complete description of this representation, and for the definition of the 
Wiener measure, to be used below. Results of this section can actually be found in |Ginj . 
but the cluster expansion used there is very intricate. 

We start with the partition function 



n^A(wi) n exp{- / U{uJi{t)-ujj{t))dty (31) 

Here, is the permutation group of elements, e = 1 for bosons and e = —1 for 
fermions, |7r| is the number of transpositions of tt, X\{u)) is 1 if u){t) G A for all ^ t ^ /? 
and is otherwise, and U{x) represents a pair interaction potential. For simplicity we 
assume that U is nonnegative and summable. Recall that the Wiener measure [Ijj 
satisfies the following property; if < < ■ • • < i^i < /3 and if / is a function W^'^ — > R, 
we have 

dT^4(^)/(u;(ti),...,u;(t„)) 

dxi . . . / dXnVti (2:1 - x)'4)t.,-t^ (X2 - Xi) . . . 'i/'/3-t„ {v - Xn)f{xi, ...,Xn), (32) 

where iptix) is the Gaussian 

V't(x) = (27rt)-'^/2 g-xV2t (33) 

The first step is to derive an expression for the partition function that is of the form 
A permutation vr G Sn can be decomposed in k cycles of lengths £i, . . .£k whose sum 
is A^. Given £i, . . . there are 

1 A^! 

corresponding permutations. Furthermore, the integrations over xi, . . . ,Xi and over wi, 
. . . , LUi, where luj is a path from Xj to Xj+i (with x^+i = xi), can be performed by a single 
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integration over x, followed by an integration over a path uj with the measure wix- We 
denote uj = {£, x, co) where £ is a positive integer, x G A, and a; is a path [0, £f3] i— > M'^ with 
a;(0) = io{£(3) = x. We consider the following measure fi 

^^+1 

df,iu) = —^dx dw^ifM xaM 

JJ expj- /'^C/(w(m/3 + t) - w(n/3 + t))dt}. (34) 

0^m<n^^-l ^ 

The partition function H31|) can be written as 
^ ^ 12 Jn [ dM(^i)-- - / dM(^iv) 

W exp|-^^/ U{ui{ml3 + t)-ujj{n(5 + t))dtY (35) 

1 C i<j < N m=0 n=0 ''^ 



We obtain (0) by setting 

C(w,^') =exp{- ^ ^ / U{uj{'ml3 + t)-u'{nl3 + t))dt^-l. (36) 

m=0 n=0 

We now establish the criterion ©. The condition below involves z, /?, and f7, see H40[) . 
and is not the most general that can be achieved. It is enough for the purpose of an 
illustration of the use of cluster expansions, however. 

We take a{u) = (— logz)^. Since 1 — e~* ^ t, we have 



e-i e'-i „^ 

U{uj{ml3 + t)-u;'{nP + t))dt. (37) 

m— n r,—n •'0 



m=0 n=0 

Now for all y and all < t' < ^'/3, it is not hard to see that 



d^' / ^w2MU{y - uj'it')) = (27r£'/3)~'^/2 f (33) 



We get then 



I dMc.'mu;,u')\z~'' ^ ^^-^ I U{x)dx X; r'^'- (39) 



The criterion is fulfilled if the right side of the expression above is smaller than (— log z)l, 
that is, if 

13 



(27r/?)'^/2 



/ U{x)dx r'^/^ ^ - logz. (40) 

This assumes in particular that z ^ 1. The thermodynamic pressure is analytic in z,/3 
in the range of parameters where (|40j) holds, and no condensation takes place, whether 
classical or of the Bose-Einstein type. 

Correlations are given by 'reduced density matrices' and require more efforts. 
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